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Abstract 

There is growing evidence that electronic and molecular networks present some common 
universal properties, among which the existence of a self-organized intermediate phase. In 
glasses, the latter is revealed by the reversibility window obtained from complex calorimetric 
measurements at the glass transition. Here we focus on amorphous networks and we show 
how this intermediate phase can be understood from a rigidity percolation analysis on size 
increasing clusters. This provides benchmarks and guidance for an electromechanical analogy 
with high temperature super conductors. 

1 Introduction 

In the recent years, more and more studies devoted to the understanding of the metal- 
insulator transition (MIT) have stressed the possible existence of a self-organized elec- 
tronic phase [0, |3[] that may explain the superconducting state. In these studies, in 
addition to the prediction or the measurement of a critical dopant concentration n c i, 
it is noteworthy to stress that although the mobile charge carrier concentration should 
tend to zero at the MIT and therefore also the interaction between the carriers, the 
temperature T c should tend to zero, in contrast with current observation. Moreover, a 
second anomalous behavior is observed at a concentration n C 2 from the superconductive 
phase to the non-superconductive Fermi liquid. This second transition is much more 
abrupt and seems first order in character Q. Illustrative examples of these behaviors 
can be found in LSCO (La2~ x Sr x CuOi) from measurements of the filling factors by 
Meissner volume (i.e. measuring the fraction of the material that is super conductive) 
and the transition temperature T c (x) with respect to Sr concentration 0. The nature 
of this intermediate phase in the dopant region n cl < n < n C 2 therefore still needs 
to be understood. However, beyond the details of modelling^, a general consensus 
emerge, which puts forward the idea of the formation of modulated spatial patterns of 
mesoscopic length scale as being responsable for HTSC|3|. ||. 

Much more simpler examples of self-organization and the existence of an inter- 
mediate phase can be provided by network glasses [0]. Here, instead of following the 
concentration of mobile carriers, focus is made on elastic deformations of an amor- 
phous network. Thus, instead of MIT, percolative stiffness transitions are obtained 
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upon changing the number of possible deformations which in turn can be related to 
the connectivity of the molecular network. In order to describe these deformations, 
one simple and elegant idea is to translate the covalent interatomic valence forces of 
the atoms in bond-stretching n" and bond-bending mechanical constraints , using 
(mean- field) Maxwell constraint counting||]. A very particular point is reached when 
the number of mechanical constraints n c = n" + nr c per atom equals the degrees of 
freedom per atom in 3D, i.e. n c = 3 (corresponding to optimal glass formation) fllCfl . 
This result has been obtained independently from rigidity percolation by Thorpe [|Tl|l 
on disordered networks, showing that the number of zero frequency solutions (floppy 
modes) / of the dynamical matrix of the network equals / = 3 — n c . The enumeration 
of constraints has been performed on a network of N atoms composed of n r atoms that 
are r-fold coordinated yielding n" = r/2 bond-stretching constraints and = 2r — 3 
bond-bending constraints for a r-fold coordinated atom. The vanishing of the number 
of floppy modes is obtained when the mean coordination number of the network reaches 



the critical value of 2.4, identified with an elastic stiffness transition) 12]. This point 



defines the transition between a floppy network which can be deformed without any 
cost in energy and a stressed rigid network which has more constraints than degrees 
of freedom. At the mean coordination number of 2.4, the network is isostatic. This 
new phase transition has been studied with considerable success in chalcogenides both 
in experiments and numerical simulations ||13||. Applications of rigidity theory have also 
been reported in various fields such as granular matter, biology and computational 
. However, in the more recent years it has been shown that there should be 
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two transitions |15j instead of the previously reported single transition suggesting that 
the mean-field constraint counting alone may be insufficient to describe completely the 
network change. Results on Raman scattering and measurements of the kinetics of the 
glass transition using modulated differential scanning calorimetry (MDSC) have indeed 
provided evidence |15|, [RJ for two vibrational thresholds thus defining a stress-free inter- 
mediate phase (isostatically rigid) in between the floppy and the stressed rigid phases, 
for which independent evidence is obtained from numerical simulations ||17|| . 

In this paper, we will construct a systematic approach that enables us to perform 
constraint counting on size increasing structures thus permitting to take into account 
medium range order effects such as small rings. In doing this, we show that these small 
rings mostly control the values of the critical coordination numbers and the width of 
the intermediate phase. The construction will also show how isostatic regions and self- 
organization influence the absolute magnitude of the width. The simplest case that 
can be build up corresponds to Group IV chalcogenide glasses of the form B x Ai_ x 
with coordination numbers ta = 2 and r# = 4 defining the mean coordination number 
f = 2 + 2x. Size-increasing cluster approximations (SICA) have been used in order to 
generate medium range order (MRO) from sets of clusters on which we have realized 
constraint counting. The results show two transitions. A first one at which the number 
of floppy modes vanishes and a second one (a "stress transition") beyond which stress 
in the structure can not be avoided anymore out of ring structures. In between, this 
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defines an almost stress-free network structure for which the rate of isostatic regions can 
be computed. We will extend these results on solid electrolytes for which percolative 
conductivity should be detected. 



2 Size increasing cluster approximations and con- 
straint counting 

Size increasing cluster approximations (SICA) have been first introduced to describe the 



medium and intermediate range order in amorphous semi-conductors |1| but their use- 
fulness has been also stressed for the description of the formation of quasi-crystals[ |i"9fl 
and fullerenes[p0|1. This approach emphasizes the rapid convergence of significant MRO 
structures to a limit value when the size of the considered clusters is increasing. Ide- 
ally, the infinite size cluster distribution would yield the exact statistics of MRO in 
the structure. The construction of these clusters is realized in Canonical Ensemble 
with particular energy levels corresponding to bond creation between short range order 
molecules (basic units) which are used as building blocks from step 1=1 (corresponding 
here to the reported mean-field approach |J) to arbitrary 1. This construction should be 



realized at the formation of the network, when T equals the Active temperature Tf [21 
Since we expect to relate the width of the intermediate phase with the ring fraction, 
we will restrict our present study to Group IV chalcogenides of the form Si x Sei_ x . For 
the latter, there is strong evidence that at the stoichiometric concentration x=0.33 a 



substantial amount of edge-sharing SiSe^ tetrahedra|22[ can be found. Therefore, 
we select basic units such as the A 2 (i.e. Se 2 ) chain fragment and the stoichiometric 
molecule (e.g. SiSe^). These basic units have respective probabilities 1—p 
and p — 2x/(l — x), x being the concentration of the Group IV atoms. We associate 
the creation of a chain-like A 2 — A 2 structure (see Table I) with an energy state of 
Ei, isostatic A 2 — BA 2 bondings with an energy gain of E 2 and corner-sharing (CS) 
and edge-sharing (ES) BA4/2 tetrahedra or any ring structure respectively with E 3 and 
E 4 . The probabilities of the different clusters have statistical weights g{Ej) which can 
be regarded as the degeneracy of the corresponding energy level and correspond to the 
number of equivalent ways a cluster can be constructed. Examples of statistical weights 
for the step I = 2 are shown in Table I and in Refs|18|, |19fl . 

Due to the initial choice of the basic units, the energy E 2 will mostly determine 
the probability of isostatic clusters since this quantity is involved in the probability of 
creating the isostatic BA± cluster (a A 2 — BA A / 2 bonding). As a consequence, if we 
have E 2 <C E\, E%, E4, the network will be mainly isostatic in the range of interest. 

The calculation of the number of bond-bending and bond-stretching mechanical 
constraints n£ per atom is performed on each cluster by Maxwell counting and redun- 
dant constraints in ring structures are removed following the procedure described by 
Thorp eflllj. It is obvious from the construction that all the cluster probabilites will 



depend only on two parameters (i.e. the factors t\je 2 and e 3 /e 2 ) and eventually e 4 /e 2 
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if one considers the possibility of edge-sharing (ES) tetrahedra or rings. One of the two 
factors has to be composition dependent since a conservation law for the concentration 
of B atoms x"> has to be fullfilled at any step / of the construction |23f : 



x (l) = x (1) 



This means that either the Active temperature Tf or the energies Ei depend on x[[21 
but here only the ei(x) dependence is relevant for our purpose. Of course, when the step 
I will increase, the number of isomers will also increase among those also the number of 
different types of rings. The construction has been realized up to the step / = 4 which 
already creates clusters of MRO size. For each step I, we have determined either ei/e2 
or ea/e2 solving equ. (|l|) and computed the total number of constraints n l c : 

Eiii N iPi 

where n c u) and Aj are respectively the number of constraints and the number of atoms of 
the cluster of size / with probability p^. Once the factors become composition dependent, 
it is possible to find for which concentration x (or which mean coordination number f) 
the system reaches optimal glass formation where the number of floppy modes / = 3— n c 
vanishes. 



3 Results 

Various possibilities can be studied within this framework. The simplest case which 
can be investigated at the very beginning is the random bonding case which is obtained 
when the cluster probabilities Pi are only given by their statistical weights (in other 
words, the factors are set to one). A single solution is obtained for the glass optimum 
point defined by / = at all SICA steps, in the mean coordination number range 
[2.231,2.275], slightly lower than the usual mean-field value of 2.4. Since there is only 
one solution, no intermediate phase in the case of random bonding is found. 

Self-organization can be obtained by starting from a floppy cluster of size I (e.g. a 
chain-like structure made of a majority of A atoms), and allowing the agglomeration of 
a new basic unit onto this cluster to generate the cluster of size I + 1 only if the creation 
of a stressed rigid region can be avoided on this new cluster. The latter occurs when 
two basic units are joined together. With this rather simple rule, upon increasing 

f one accumulates isostatic rigid regions on the size increasing clusters because BA^/ 2 
units are only accepted in A 2 — BA4/2 isostatic bondings with energy E 2 . Alternatively, 
we can start from a stressed rigid cluster which exist at higher mean coordination 
number (f < 2.67) and follow the same procedure but in opposite way, i.e. we allow only 
bondings which lead to isostatic rigid regions, excluding systematically the possibility 
of floppy A 2 — A 2 bondings. In the case of self-organized clusters, the simplest case to 
be studied is the case of dendritic clusters, when ring creation is avoided. For an infinite 
size I, this would permit to recover the results from Bethe lattice solutions or Random 
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Bond Models[£4] for which rings are also exluded in the thermodynamic limit [£5|]. A 
single transition for even / steps at exactly the mean-field value f = 2.4 is obtained 
whereas for the step 1 = 3, there is a sharp intermediate phase defined by / = (still 
at f = 2.4) and the vanishing of floppy regions (i.e. ei/e 2 is zero) at f = 2.382(6). 
The probability of floppy isostatic rigid and stressed rigid clusters as a function of the 
mean coordination number has been computed and shows that the network is entirely 
isostatic at the point where / = (solid line, fig. 1). 

The intermediate phase shows up if a certain amount of medium range order (MRO) 
is allowed. This is realized in the SICA construction by setting the quantity e^j e 2 non 
zero, i.e. edge-sharing (ES) tetrahedra BA^/2 leading to four-membered rings B 2 A± can 
now be created at the growing cluster steps. Two transitions are now obtained for every 
SICA step. The first one is still located at f c \ = 2.4 (see fig. 1). There, the number 
of floppy modes / vanishes. The second one is located at f c2 . When starting from a 
floppy network close to f = 2 and requiring self-organization (in this case, allowing 
only floppy or isostatic bondings), this point corresponds to the network composition 
beyond which stressed rigid regions outside of ring structures, which are created by the 
dendritic connection of at least two BA4/2 units can not be avoided anymore. We call 
this point the "stress transition". We show in fig. 1 the I = 2 result where / = at 
f cl = 2.4 for different fractions of ES tetrahedra, defining an intermediate phase Af . It 
is noteworthy to stress that the first transition at f c \ corresponding to Phillips' glass 
optimum does not depend on the ES fraction, as well as the fraction of stressed rigid 
clusters in the structure. To ensure continuous deformation of the network when B 
atoms are added and keeping the sum of the probability of floppy, isostatic rigid and 
stressed rigid clusters equal to one, the probability of isostatic rigid clusters connects 
the isostatic solid line at f c2 . Stressed rigid rings first appear in the region T c \ < f < f c2 
while chain-like stressed clusters (whose probability is proportional to e 3 ) occur only 
beyond the stress transition, when e^, ^ 0. This means that within this approach, when 
f is increased, stressed rigidity nucleates through the network starting from rings, as 
ES tetrahedra. It appears from fig. 1 that the width Af = f c2 — f c \ of the intermediate 
phase increases with the fraction of ES. We have represented the width as a function 
of the overall MRO fraction at the rigidity transition in fig. 2 which shows that Af 
is almost an increasing function of the ES fraction as seen from the result at SICA 
step I = 4. Here, there is only a small difference between allowing only four-membered 
rings (ES) (lower dotted line) or rings of all sizes (upper dotted line) in the clusters. 
Concerning the nature of the structure in the intermediate phase, one can see from fig. 
1 and the insert of fig. 2 that the probability of isostatic clusters is maximum in the 
window Af , and almost equal to 1 for the even SICA steps, providing evidence that the 
molecular structure in the window is almost stress-free. Furthermore, the latter can be 
extremely broad if a high amount of edge-sharing is allowed. For a fraction of ES 77 
is close to one, molecular stripes emerge in the structure which are isostatic and every 
addition of new B atoms lowers the amount of remaining A-chains by converting them 
in ES BA4/2 stripes. These molecular stripes are stress-free only if their size is infinite 
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11 1, which is the case for rj ~ 1 in the thermodynamic limit. 



4 Application to chalcogenides 

As mentioned above, chalcogenide glasses are the first systems which have been carefully 
studied in the context of self-organization. SICA provides therefore a benchmark to 
check the results obtained. To be specific, Raman scattering has been used [^7[ ^8 



as a probe to detect elastic thresholds in binary Si x Se\^ x and Ge x Sei_ x or ternary 



Ge x As x Se\_2x glasses p9| . Changes in the germanium or silicon corner sharing mode 
chain frequencies have been studied with mean coordination number f of the glass 
network. These frequencies exhibit a change in slope at the mean coordination number 
f cl = 2.4 and a first order jump at the second transition f c2 . In germanium systems, 
the second transition is located around the mean coordination number of 2.52 whereas 
f C 2 = 2.54 in Si based systems. For both systems, a power-law behavior in f — f C 2 
is detected for f > f c2 and the corresponding measured exponent is very close to the 
one obtained in numerical simulations of stressed rigid networks [2g] . Moereover, a clear 



correlation between these results and the vanishing of the non-reversing heat flow AH nr 
(the part of the heat flow which is thermal history sensitive) in MDSC measurements 
has been shown p8| . 

The SICA approach shows that the width Af of the intermediate phase increases 
mostly with the fraction of ES tetrahedra. We stress that the width should converge 
to a lower limit value of Af compared to the step 1 = 2, therefore one can observe the 
shift downwards when increasing / from 2 to 4. This limit value is in principle attained 
for I — > oo, or at least for much larger steps than / = 4 fl30 |. For Si-Se, Af = 0.14 is 



somewhat larger than for Ge-Se (Ar = 0.12) consistently with the fact that the number 
of ES is higher in the former |28j . 

5 Application to solid electrolytes 

One interesting field of application of cluster construction and constraint counting algo- 
rithms is the area of fast ionic conductors (FIC) |3l|, |32| , which has recently gained at- 
tention because of potential applications of these solid electrolytes in all solid state elec- 
trochemical devices and/or miniaturized systems such as solid state batteries. This has 
been made possible by substitution of more polarizable sulfur atoms, replacing the oxy- 
gen in usual oxide glasses and has led to a substantial increase of the conductivity [[33 



as 



high as 10 3 Q 1 .cm 1 , three orders of magnitude higher than the conductivity of anal- 
ogous oxide glasses EM. Surpisingly, the extension of constraint theory from network 



glasses (as Si x Sei_ x ) to FIC has been only reported for a few oxide glasses|35], [3(| . 
As a consequence, percolative effects have not been studied so far with the network 
change in FIC, although it seems fundamental for the understanding of the mobile al- 
kali cations' motion. Since the conductivity of the semi-conductor is proportional to 
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the mobility fi, which in turn is related to the deformation of the network ||37||, and the 
free carrier concentration ni by o = /itilg, one may expect that the mobility in a floppy 
FIC should be substantially higher compared to the cation mobility in a stressed rigid 
network. This means that the conductivity a should display some particular behavior 
in the stress-free intermediate phase and at the two transitions. 

One can extend the SICA approach to the present case by considering the simplest 
binary conducting glass, which is of the form (1 — x)SiX2 — xM 2 X with X an atom of 
Group VI (X=0,S,Se) and M an alkali cation (M=Li,Na,K,. . . ). The local structure 
can be obtained from NMR measurements and is mainly made of so-called Q 3 and Q 4 
units |38|. The former corresponds to the usual silica tetrahedron made of one silicon 
and four Group VI atoms at the corner [e.g. SiSe^/ 2 ] while the latter has one oxygen 
atom ionically bonded to the alkali cation [e.g. SiSef, 2 Na®] (fig. 3). 

Then, the probabilities can be evaluated for different steps of cluster sizes following 
the procedure described in the previous sections. It appears that the creation of a 
Q 4 — Q 4 connection leads to a stressed rigid cluster if the number of constraints is 
computed, while the Q 4 — Q 3 and Q 3 — Q 3 connections yield respectively isostatically 
stressed and floppy clusters. The SICA results show again that the intermediate phase 
shows up if a non-zero fraction of small rings is allowed in the structure, as displayed 
in fig. 4. Here in contrast with chalcogenides, the glass optimum correponding to the 
vanishing of the number of floppy modes respresents the upper limit of the intermediate 
phase, which is consistent with the fact that one starts from an almost stressed rigid 
molecular network at low modifier concentration. Increase of the alkali content leads to 
an increase of Soppiness. In the oxide system (1 — x)Si02 — xM 2 0, the width Ax should 



be very small or zero since the fraction of ES in the oxide systems is almost zero[21 



Still, percolative effects are expected at the concentration x = 0.2 corresponding to 
the transition from rigid to floppy networks. In sulfur and selenide glasses such as 
(1 — x)SiS 2 — xNa 2 S, the width should be much broader because of the existence 



of a high amount of edge-sharing tetrahedra in the SiS 2 or SiSe 2 base networks |39 
In the sulfur base glass, 29 Si NMR have shown that the rate of ES should be about 



0.5, slightly higher than in the selenide analogous system ||(|. As a result, one should 
observe a window of about Ax = 0.09. Unfortunately, results on these systems are only 
available for an alkali concentration x > 0.2f4T|. However, in the different silica based 
glasses, a rigidity transition has been observed |[42|| at the concentration x = 0.2 which 
should provide guidance for forthcoming studies in this area. 

Finally, temperature effects should be observable close to this transition. Since 
the concentration of alkali free carriers depends on the temperature (the higher 
the temperature, the higher til), an increase of the temperature T should lead to a 
decrease of the number of network constraints, the fraction of intact bond-stretching 
constraints n™ of the alkali atom being proportional to 1 — n^. Consequently, a shift 
of the mechnical threshold (/ = 0) to the higher concentrations should result from an 
increase of T. 
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6 Summary 



In this article, we have shown how stress change in molecular systems can be understood 
from the combination of cluster construction and constraint counting. We have found 
that there is a single transition from floppy to rigid networks in a certain number of 
structural possibilities. An intermediate phase intervenes when a fraction of non-zero 
molecular isostatic stripes are allowed, which makes this new phase globally stress-free. 
The similarity of this new phase with the superconducting filamentary phase in HTSC[H] 



but also in biological networks JLJ] is striking and should be worked out in more depth 
in the future. 
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symbol 


mechanical nature 




probability 


n c 


A 4 


floppy 




4(l-p) 2 ei 


2.0 


BA 4 


isostatic 


o 


16p(l - p)e 2 


3.0 


B 2 A 4 


stressed 




16p 2 e 3 


3.67 


B 2 A 4 


stressed 




72p 2 e 4 


3.33 



Table 1: The four clusters generated at step I = 2 with their unrenormalized probabil- 
ities and their number of mechanical constraints. The factors are Boltzmann factors 
involving bond energies and the Active temperature, et = exp[— Ei/Tj] 
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T> 0.6 
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O 
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stressed 



floppy \. 



Ar 



_l 



2.2 2.4 

Mean coordination number r 



2.6 



Figure 1 : Probability of finding floppy, isostatic rigid and stressed rigid clusters as a function 
of mean coordination number f at SICA step I = 2 for different fractions of edge-sharing. The 
solid line corresponds to the dendritic case where no rings are allowed (e^ = 0). The broken 
lines correspond to situations whith a fraction of ES at the stress transition of ??=0.156, 0.290 
and 0.818. The filled squares indicate the point f C 2 at which the stress transition occurs in 
the case r] = 0.156. 
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Fraction of MRO at the rigidity transition 

Figure 2: Width of the intermediate phase Af as a function of the fraction of ring clusters 
at the rigidity transition for / = 2 (solid line), I = 3 (dashed line) and I = 4 (dotted lines). At 
step I = 2, the non-zero width comes only from the edge-sharing BA±/ 2 tetrahedra. For / > 2, 
different rings sizes (4, 6, 8) have been taken into account. The lower dotted line corresponds 
to the intermediate phase at I = 4 if only ES are allowed as stressed rigid fragments. The 
upper dotted line is the same quantity but allowing any ring structure. The insert shows the 
probability of isostatic clusters with mean coordination number f for I = 4 (dotted line) and 
I = 3 (dashed line). The shaded region of I = 4 represents the intermediate phase. 




Figure 3: The two local structures of the fast ionic conducting (1 — x)SiX2 — XM2X, {X = 
O, S, Se; M = Li, Na, K) glasses, denoted as Q A and Q 3 units, with respective probabilities 
1 — p and p = 2x/(l — x). 
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Figure 4: Width of the intermediate phase Ax as a function of the rate of edge-sharing (ES) 
in the base glass. The insert shows the solutions x c \ and x C 2 obtained from SICA analyis, 
again as a function of ES. 
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